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The dipole trap model is able to explain the main properties of the apparent metal-to-insulator 
transition in gated high mobility Si-inversion layers. Our numerical calculations are compared with 
previous analytical ones and the assumptions of the model are discussed carefully. In general we 
find a similar behavior but include further details in the calculation. The calculated strong density 
dependence of the resistivity is not yet in full agreement with the experiment. 
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I. INTRODUCTION 

Since it's discovery in 1995 1 , the metal-insulator tran- 
sition in two dimensions (2D) was investigated carefully 2 , 
as it's finding is in apparent contradiction to the scaling 
theory of localization 3 . According to the latter, in the 
limit of zero temperature, a metallic state exists only in 
three dimensions, but in two dimensions disorder should 
always be strong enough in order to lead to an insulating 
stated. Nevertheless, high-mobility n-type silicon inver- 
sion layers showed for high electron densities a strong 
decrease of resistivity p towards temperatures below a 
few Kelvin, manifesting the metallic region, and a strong 
exponential increase of the p for low densities demon- 
strating the insulting regime. A very similar behavior 
was observed in many other semiconducting material sys- 
tems at low temperatures. 

After the unexpected finding, several models were 
suggested in order to explain the metallic behavior 
in 2D. The most important ones are i) temperature- 
dependent screening^£, ii) quantum corrections in the 
diffusive regim o 7 i 8 i 9 , iii) quantum corrections in the bal- 
listic regim o 10 : 11 , and iv) scattering of electrons according 
to the dipole trap modeli^. As there are argumentations 
for all that different models in the literature, we do not 
want to repeat them here in detail. A clear decision for 
one of the suggestions could not been drawn yet and fur- 
ther work on the models has to be carried out. 

The dipole trap model was introduced by Altshulcr 
and Maskr\*i£ (AM) especially for Si-MOS structures, as 
it is known that the misfit at the silicon/silicon-oxide 
interface produces charged defect states in the thermally 
grown oxide layeoiiiiiii^. AM could show that a hole 
trap level at energy E t which is either filled or empty, 
depending on it's position relative to the Fermi energy 
Ep, can lead to a critical behavior in electron scattering 
if E t and Ep are degenerated This dipole trap model 
can explain the main properties of the metal-insulator 
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transitions in gated Si-MOS structures^. 

In this work, we present numerical calculations of the 
temperature and density dependent resistivity due to 
electronic scattering in the dipole trap model. With these 
calculations we are able to checking the analytical calcu- 
lations with it's approximations. Due to the numerical 
procedure, we can include further details and and inves- 
tigate their influence. 

For the analytical calculations AM made a number of 
assumptions. These are: al) the trap states possess a 5- 
like distribution in energy, a2) the spatial distribution in 
the oxide is homogeneous, a3) the occupied states behave 
neutral and cause no scattering of 2D electrons whereas 
the unoccupied states are positively charged and lead to 
scattering (hole trap), a4) a charged trap state is effec- 
tively screened by the 2D electrons so that the result- 
ing electrostatic potential can be described by the trap 
charge and an apparent mirror charge with opposite sign 
on the other side of the interface, a5) the scattering effi- 
ciency of the 2D electrons is described by a dipole field 
of the trap charge and it's mirror charge, a6) a parabolic 
saddle point approximation for the total potential of the 
trap states was used in order to perform analytical cal- 
culations, a7) the energy of the trap state E t Q is fixed 
relative to the quantization energy Eo of the 2D ground 
state inside the inversion potential, and a8) the Fermi 
energy Ep in the 2D layer is either independent of or is 
the same as in the 3D substrate. 

In contrast to AM, our calculations were performed nu- 
merically, so that several limitations of their calculations 
could be dropped. Our improvements concern il) the de- 
tailed spacial dependence of the electrostatic potential is 
taken into account instead of the parabolic saddle point 
approximation, i2) the energy of the trap state E t o is 
fixed relative to the conduction band edge Eq b ■ As a re- 
sult of our calculations, we find a similar behavior of the 
calculated resistivity as AM and we calculate in addition 
the density dependence of the resistivity. 

According to the restricted space in the original AM 
work, some of the used equations were not derived there. 
We will discuss these equation and considerations in de- 
tail in the main part. For better readability of our paper, 
some details were put into appendices. Please note that 
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we will use SI units throughout this work. 



II. MODEL CONSIDERATIONS AND 
NUMERICAL CALCULATIONS 

The misfit at the Si/SiC>2 interface layer leads to dif- 
ferent kinds of defects and trap statesi^^i^^. In the 
considered AM model it is assumed that a relative large 
number of hole trap states exists. If such a trap state 
captures a hole, it is positively charged, otherwise it is 
neutral. In Fig.EJthe trap state is depicted schematically 
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The last term in Eq. pleads to a down bending of the en- 
ergetic position towards the interface and causes a max- 
imum in the total trap energy E t (z) as shown in Fig. 
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FIG. 1: Schematic representation of the trap state together 
with the 2D electron system (2DES). For E t > Ef the trap 
state is positively charged and scatters electrons in the 2D 
layer whereas for Et < Ef the trap is neutral, i.e. it is occu- 
pied with an electron and does not act as a scattering center. 
Note that the main recharging effect comes from the ener- 
getic position E t of the trap state. It varies strongly with the 
applied gate voltage V g , whereas the Fermi-energy Ef shows 
only small variations with changing electron density on the 
same scale. 



As described in the introduction, it is further as- 
sumed that al) all trap states exist at the same en- 
ergy Eto if no external field is applied and posses a2) 
a spatially homogeneous density distribution in the ox- 
ide layer. But a potential gradient due to an applied gate 
voltage V g causes a linear increase of the trap energy po- 
sition Et = Eto + eVgz/d, where z is the distance from the 
Si/SiC>2 interface (z < 0) and d is the distance between 
gate electrode and that interface (i.e. the thickness of the 
oxide layer). 

For the electrostatic potential inside the oxide layer, 
also the screening effects of the inversion layer have to be 
taken into account. For 2D electrons in a Si-(OOl) layer, 
the screening radius is equal to ag/4. If the trap dis- 
tance from the interface \z\ exceeds the screening radius, 
the in-plane components of the electrostatic field caused 
by the charged trap will effectively be screened. In that 
case, the electric field and the potential in the oxide can 
be described by the trap charge and an apparent mirror 
charge with opposite sign on the other side of the inter- 
face (assumption a4) . The potential of the charged state 
caused by the image charge is $ = e/(2 • 4ireoe ox z) in SI 
units with e = 8.854 x 1CT 12 Fm _1 and e ox w 3.9, the 
relative dielectric constant of the oxide. Thus the total 
energy of the charged trap state can be given as 



E t {z) = E, 



to 



8ir£ S ox Z 



for z < . (1) 



FIG. 2: Schematic representation of energies in the dipole 
trap model. The dashed blue line on the left represents the 
trap energy without, the full blue line with mirror charge 
potential. Note that the distance scale on the left and right 
hand side of the interface is chosen differently in order to 
increase visibility. 



The trap charge together with it's mirror charge form 
a dipole perpendicular to the interface plane. Thus, for 
distances larger than 2\z\, the scattering potential expe- 
rienced by the 2D electrons can be described by a dipole 
field which falls of with 1/r 3 (assumption a5). This is 
in agreement with the long range field of a screened 
Coulomb potential in two dimensions and leads to a 
consistent description. AM have calculated the clas- 
sical scattering cross section for momentum relaxation 
a m (E, z) for such a dipole field for electrons with kinetic 
energy E as 



X (E, z) = 2.74(e 2 z 2 / 



iire e*E) 



1/3 



(2) 



with e* = (e ox + £si)/2 ~ 7.9, the effective dielectric 
constant for the 2D electron system (2DES). 

Whether a trap state is charged or not, depends on it's 
energetical position relative to the Fermi energy Ep (as- 
suming thermal equilibrium for the occupation) . The oc- 
cupation function corresponds to a modified Fermi-Dirac 
distribution, where the degeneracy of empty and filled 
states is taken into account. AM have assumed that the 
(positively) charged trap state can have either spin up or 
down and is thus doubly degenerate, while the neutral 
state has no degree of freedom and is not degenerate. 
From that the probability of a trap state to be charged 
follows as 



p+(z) 



1 



■ exp 



(3) 



with k the Boltzmann constant. 
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For the occupation of the trap states only the relative 
position of the trap energy E t (z) to the Fermi energy Ep 
is important. But the difference Ep — E t (z) can not be 
derived directly - it has to be calculated from the two 
individual energies which depend on different variables. 
According to Eq. ^ the z-dependence of the trap energy 
can be calculated, but one has to fix it's zero-position 
E t Q. AM have assumed (a7) a fixed energetical distance 
of the trap state relative to the quantization energy E 
of the electronic ground state in the (nearly triangular) 
inversion potential. But E Q depends on the strength and 
shape of the inversion potential and via electron-electron 
interaction on the 2D electron density n s . Thus it seems 
not realistic that the energy of the trap state is fixed rel- 
ative to Eq, but rather that it is fixed relative to the en- 
ergetic position of the conduction band edge Eq b (which 
is our improvement i2). 

Equation ^ can be used as given, by noting that the 
energy is defined relative to the conduction band edge 
Ecb- On the other hand the ground state energy E 
has to be calculated for the inversion potential, which 
itself depends on n s and the depletion charge N~ and by 
including the electron-electron interaction^. As Ep — Eo 
follows from the electron density together with Eq 
one gets the position of Ep relative to the conduction 
band edge and the difference Ep — E t (z) can be used for 
p+(z) in Eq.El 

In the Drude-Boltzmann approximation, the electrical 
resistivity p, equal to the inverse conductivity a 



follows by calculating the effective transport scattering 
time f. The detailed calculation is performed in Ap- 
pendix [S] 

As a result one gets that 

~ = N+v(E)a m {E,z m ), (5) 

T 

can be expressed by the effective values Nj~ s for the num- 
ber of charged trap states per area, v(E) the electron 
velocity, a m (i?,z m ) the scattering cross section, and the 
average electron energy E as given in Appendix^ These 
effective values depend on all the important variables of 
the systems, i.e. on T, V g and so on. By inserting Eq.[S] 
into Eq.QJ one gets already the dependence of the resis- 
tivity p on the different parameters 

P = ^ N+ v(E) a m {E, z m ) , (6) 

and we have verified equation Eq. 7 in Ref. [T3 . 

From here our treatment of the subject is quite dif- 
ferent from that of AMi£. They have evaluated Eq.|S| 
analytically whereas we perform the calculation of it nu- 
merically. But in order to be able to solve Eq.^J AM 
have used a parabolic (saddle-point) approximation for 
the z-dependence of the trap energy (assumption a6). 



The analytical expression of AM (Eq. 9) contains a tem- 
perature independent prefactor pa and a temperature de- 
pendent scaling function R(V g ,T). For their case (A) of 
the temperature dependence of Ep, they get a critical 
behavior with R(T) increasing for V„ > Vg and decreas- 
ing for V g < Vg (see Fig. 1 in Ref. . similar to what 
is observed experimentally. For their case (B), R(T) is 
always increasing with temperature and no critical be- 
havior comes out. 

In our numerical evaluation of Ea. lA6l we use the exact 
dependence of E t (z) as given by Eq.n(our improvement 
il). The numerical treatment prevents errors due to the 
parabolic approximation, but enables us also to include 
further details, which also cannot be solved analytically. 

AM have calculated the temperature dependence of 
the resistivity in close vicinity of the critical density 
where the behavior changes from metallic to insulating 
behavior. We have calculated also the direct density de- 
pendence over a larger range for different temperatures. 
Fig-El shows how the effective number of charge trap 
states N^g depends on n s . 




FIG. 3: Effective number of charged trap states per area iVi 
vs. electrons density (i.e. varying gate voltage) at tempera- 
tures of T = 1, 3, and 5 K. 



As can be seen there is a very strong variation above 
n s pa 1.5 x 10 11 cm" 2 , where the maximum of the trap en- 
ergy Et(z) is just degenerate with Ef. This strong varia- 
tion comes from the fact that as soon as the maximum of 
E t (z) is below Ep only an exponentially small number of 
traps is still excited (i.e. charged) and the scattering effi- 
ciency decreases accordingly. As p is nearly proportional 
to ATg, such strong variations have not been observed 
experimentally. This discrepancy to the experiment can 
possibly be explained that in real 2D Si-MOS structures 
either the trap states do not have a <5-like distribution in 
energy or that in addition other scattering sources exist. 
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III. CONCLUSIONS 

We have shown that the numerical calculations of the 
temperature dependent resistivity give similar results as 
the analytical methods by AM. The strong density de- 
pendence of ATg and thus of p which follows from the 
calculation is not in agreement with experimental find- 
ings. In order to possibly resolve this discrepancy further 
calculations should be performed within the dipole trap 
model. The numerical procedure allows incorporation of 
further effects and realistic assumption like energetical 
broadening of the trap level, special spatial distributions 
of the defects, and detailed screening dependence. 
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APPENDIX A: TRANSPORT EQUATIONS 

The effective transport scattering time f in the Drude- 
Boltzmann approximation follows from^AS 



J dET(E)Edf/dE 
JdEEdf/dE 



(Al) 



with t(E) being the energy dependent scattering time 
and df jdE the first derivative of the Fermi-Dirac distri- 
bution function /. 

The transport scattering time t(E) has to be calcu- 
lated by integration over the individual scattering rates 



(both being three dimensional volume densities) and 
v(E) the electron velocity. 

By inserting the all expressions into Eq. IA2I one gets 



l/r(E) = c'N+zTE^ 



(A3) 



with the prefactor c' = 2.74(e 2 /8ire e*) 1/3 ^2/m*, an 
effective number of positive trap states per area = 
N t 3 (Az) eS , the effective width of positive charge layer 
(Az) eff = J dzp + (z)(z/z m ) 2 / 3 and the position z m = 
y/ed/8neQe ox Vg of the energetical maximum of the trap 
energy. 

By inserting Eq. IA3I into Eq. IA1I one gets 



/>oo 

f oc 1/E F / dEE 5/e df/dE. 
Jo 



(A4) 



Further an effective energy E can be defined so that 
formally Ea. lA3l can be preserved for the effective f, i.e. 
1/f = c' N^ffZm 3 E 1 / 6 . A simple calculation gives 



E = Ei 



5/6 



df/dE 



(A5) 



By further replacing the first derivative of the Fermi- 
Dirac function by the identity df/dE = — f(l — /) = 
-(4fcTcosh 2 ((.E - E F )/2kT))- 1 one obtains the same 

M 



expression as Eq. 8 in Ref. 

With these relations, the resistivity can exactly be 
written in terms of the effective energy E as 



P = — j N eff v ( E ) a m(E, z m ) , 



ne 



(A6) 



which corresponds to Eq. 7 in Ref. [lj, but the individual 
terms are rewritten according to our definitions above. 
A comparison with Eq.0] gives exactly 



1/t(E)= dzN+(z)v(E)a m (E,z) (A2) 



with the density of charged traps N^(z) = N t 3P+(z), 
N t3 the spatially constant density of existing trap states 



^=N+v{E)a m (E,z m ), 

T 



(A7) 



and shows that Ea. lA2l can also be rewritten for effective 
values. 
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